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The Very Massive and Hot LMC Star VFTS 682: Progenitor of a 

Future Dark Gamma- Ray Burst? 
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O ■ ABSTRACT 
(N" 

O-i VFTS 682, a very massive and very hot Wolf-Rayet (WR) star recently dis- 

covered in the Large Magellanic Cloud near the famous star cluster R136, might 
be providing us with a glimpse of a missing link in our understanding of Long 
Gamma- Ray Bursts (LGRBs), including dark GRBs. It is likely its properties 
^ i result from chemically homogeneous evolution (CHE), believed to be a key pro- 

^ ! cess for a massive star to become a GRB. It is also heavily obscured by dust 

i extinction, which could make it a dark GRB upon explosion. Using Spitzer data 

Jp"! we investigate the properties of interstellar dust in the vicinity of R136, and ar- 

O I gue that its high obscuration is not unusual for its environment and that it could 

^ i indeed be a slow runaway ("walkaway") from R136. Unfortunately, based on its 

current mass loss rate, VFTS 682 is unlikely to become a GRB, because it will 
^ . lose too much angular momentum at its death. If it were to become a GRB, it 

>■ I probably would also not be dark, either escaping or destroying its surrounding 

^ " dusty region. Nevertheless, it is a very interesting star, deserving further studies, 

O I and being one of only three presently identified WR stars (two others in the Small 

^ I Magellanic Cloud) that seems to be undergoing CHE. 

Subject headings: stars: mass-loss — stars: Wolf-Rayet — gamma-ray burst: 
general — Magellanic Clouds 

1. Introduction 

Ideally, in order to fully underst and the mapping betw een massive star progenitors and 



when and how they explode (or not, iKochanek et al.ll2008l ). we would have extensive multi 



wavelength data obtained for many such explosions both BEFORE and AFTER the event. 
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The AFTER part has certainly undergone an explosive growth in the last decade or so, due 
to many successful supernova (SN) searches such as the galaxy-targeted Lick Observatory 
Supernova Search (LOSS), the Catalina Real-time Transient Survey (CRTS), the Robotic 
Optical Transient Search Experiment (ROTSE), the Palomar Transient Factory (PTE); and 
gamma-ray burst searches such as the High Energy Transient Explorer-2 (HETE-2) and 
Swift. The BEEORE part is naturally limited by the fact that massive stars, while relatively 
bright, are only observ able in the local Un iverse {d < 10 Mpc). Here we have had to rely on 



proxi mity (SN 1987A, iKunkel et al.l 119871 ) . luck based on archival data (e.g^ 



2009 ). or the first systematic campaign to monitor future SN progenitors (ISzczygiel et al. 



Smartt et al. 



201ll ). where many nearby galaxies are observed to sufficient depth so eventually they will 



provide progenitor data for a significant number of future SNe. We have by now compiled a 



1987A ( 


White k MalinI 


198 


7) or 


(Prieto et al. 


20081 


Prieto 


2008) 



However, for very rare explosive events, such as long-duration gamma-ray bursts (LGRBs) 
the systematic approach that works for the normal core-collapse SNe is not yet possible be- 
cause the events are so rare. Therefore the prospect of future systematic observational efforts 
that could identify nearby LGRB progenitors is very dim. Here we have to rely on a more 
extended chain of reasoning, with data being supplemented with reasonable, theoretical 
guesses. Discovery of the connectioii between LGRBs and broad-lined Type Ic SNe (e.g., 
Stanek et al.ll2003l : iHjorth et al.ll2003l ). thought to result from core-collapse of hydrogen- free 
massive stars, favors two p ossible progenitor models: si ngle massive Wolf-Rayet (WR) stars 

9851). or lower mass helium stars stripped 



with rapidly rotating cores ( Eilippenko fc SargentI 



by a close binary companion (jPodsiadlowski et al. 



20041 ) . The prompt emission of an initial 



LGRB, lasting seconds or minutes, is usually followed by a multi- wavelength afterglow which 
lasts days to even years. Although X-ray afterglows of LGRBs are nearly always detected by 
Swift and Fermi, detection of optical and infrared afterglows is less common. A dark GRB is 
defined by either an absent or faint optical afterglow relative to its X-ray emission. The rapid 
detection of X -ray afterglows with Swift revealed that the dark fraction of LGRBs is about 
30% of GRBs ( lAkerlof fc Swanll2007l : iPerley et al.ll2009[ ). In general, besides the intrinsically 
optically faint GRBs, the optical attenuation of dark GRBs can be caused by dust extinction 
in GRB host gal axies, foreground e xtinction, or Lyman-a absorption by neutral hydrogen 
at high redshifts ( jPerley et al.l 120091 ). 



Study of massive stars in very nearby galaxies can supply the other missing links in 
our understanding of LGRBs. The recently uncovered very massive stars up to 3OOM0 in 
the Large Magellanic Cloud ( LMC) young cluster R 136 extend our knowledge of massive 
star formation and evolution (ICrowther et al.ll2010[ ). In particular, the newly discovered 
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very massive WR star VFTS 682, located 30 pc away from R136, drew our attention, mainly 
because of its very high fo reground dust extinction and possibility of being a GRB progenitor 
( iBestenlehner et al.ll201ll ). Its unusually high effective temperature can be understood as the 
consequence of chemically homogeneous e volution (CHE), which was proposed as the crucial 
part of the process of producing LGRBs (lYoon fc Langer II2005I ). If VFTS 682 will indeed 
make a GRB at the end of its evolution, the high foreground dust extinction could mean that 
it will become a dark GRB. However, its high mass and strong mass- loss seems to prevent 
it from being a GRB progenitor, and the potential afterglow might destroy the surrounding 
dust even if VFTS 682 eventually produces a GRB. Therefore, it is worthwhile to look into 
the details of VFTS 682's evolution and the fate of its circumstellar environment. Section 
|2] describes the dust properties in the VFTS 682 vicinity region. Section |3] investigates the 
possibility of VFTS 682 producing a LGRB and Section S] discusses the fate of dusty clouds 
around VFTS 682. Conclusions are presented in Section 



Investigating the Properties of Interstellar Dust in the Vicinity of VFTS 

682 



Among the many striking properties of VFTS 682 discussed in IBestenlehner et al. 



( 120111 ). the high value they derived for the foreground extinction to that star. Ay = 4.45 ± 
0.12, particularly drew our attention. Indeed, for a random line of sight to an object in 
the LMC, a typical value o f the interstellar extinction would be much lower. Ay ~ 0.3 



Pejcha fc Stanekl 120091 ). At the same time, other very massive, nearby stars in R136- 



(e.g. 

the central cluster of the 30 Doradus region — ^^ have estimated interstellar extinctions near 



Av ~ 2.0 (see Table 3 in Crowther et al. 



20101 ) ■ significantly lower than VFTS 682. This 



high value of extinction could provide an interesting clue to the origin of VFTS 682, because 
if such values were indeed very rare in the 30 Dor region, it would argue against VFTS 682 
being a slow runaway ("walkaway") from the R136 cluster. 

To investigate the properties of interstellar dust extinction in the vicinity of VFTS 682, 
we analyzed the archival Spitzer Space Telescope IRAC data for the 30 Dor region, which were 
obtained in 2003 by the "Comparative Study of Galactic and Extragalactic HII Regions" 
program (PI: Houck; Program ID: 63). For the purpose of this investigation, we limited our 
analyses to a region 20' x 20' centered on R136. 

In Figured] we show 3.6, 4.5, 5.8 and 8.0 /im IRAC images of the region near VFTS 682 
and 30 Doradus. We see a progression of still seeing significant stellar light in the 3.6/im 
band, so the cluster still features prominently in this band, to being completely dominated 
by warm dust emission and PAH emission in the 8.0/im band, where the cluster virtually 
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Fig. 1. — Spitzer-IRAC archival images of the 30 Doradus region in 3.6, 4.5, 5.8 and 8.0 fim 
(upper left to bottom right) IRAC photometric bands. The location of VFTS 682 is marked 
with the small circle, while the location of the R136 cluster is shown with the large circle. 
The 60" scale bar corresponds to roughly 14.5 pc for an assumed LMC distance of 50kpc. 
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30 Dor: Spitzer data 
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Fig. 2. — Mid-IR color-magnitude diagrams for the vicinity of VFTS 682 obtained using 
Daophot on the archival Spitzer IRAC images of [4.5], [3.6] -[4.5] CMD (left), and [4.5], [4.5]- 
[5.8] CMD (right). The small dots show all matched point sources in the 20' x 20' field, while 
the bigger open circles represent sources within 3' (roughly 44 pc) of R136. The red filled dot 
represents our photometry for VFTS 682, while the open star shows the values as measured 
by the SAGE project (Meixner et al. 2006). Also shown in the left panel are two reddening 
vectors corresponding to Ak value of 1.0 (about twice as high as derived by Bestenlehner et 
al. 2011 for VFTS 682), for two different values of the interstellar reddening law, Ry = 3.1 
and 4.7. For discussion see Section |2l 
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disappears. If we can use the 8.0/im band emission as a proxy for dust column density, then 
the high value of Ay towards VFTS 682 is not at all unusual in the 30 Dor region, and 
there should be lines of sight with 5 to 10 times higher values of interstellar extinction near 
the location of VFTS 682. Therefore, despite our initial expectation, the high value of Ay 
does not bring any new information about the nature of VFTS 682, namely was it ejected 
from R136 or was it born in situ. However, the presence of such large amounts of spatially 
complex dust raises the interesting possibility of there being other stars like VFTS 682 in 
the vicinity of R136, that are hidden behind still more dust. 

For a more quantitative analyses of the IRAC data, we used Daophot Jstetsonlll992l ) to 
identify point sources and measure their fluxes in the 3.6, 4.5 and 5.8 /im images. In Figure [2] 
we show the resulting mid-IR color-magnitude diagrams (CMDs). In the [4.5], [3.6] — [4.5] 
CMD we more or less agree with SAGE that VFTS 682 has a small 4.5/im flux excess, but 
we find a bluer [4.5] — [5.8] color, most likely due to differing treatments of the extended 
dust emission near VFTS 682. In any case, the amount of resolved mid-IR emission at the 
position of VFTS 682 is very small compared to its total bolometric luminosity. It is clear 
from these CMDs that the mid-IR properties of VFTS 682 are not at all unusual compared 
to other stars in the 30 Dor region, and there are many other stars nearby that are either 
more obscured or have significantly more mid-IR emission around them. 



3. Can VFTS 682 be a GRB Progenitor? 



VFTS 682 has a high temperature T^s = 52.2 ± 2.5 kK and luminosity log(L/L0) = 
6.5 ±0.2, placing it blueward of the zero-age main sequence in the Hertzsprung- Russell (HR) 
diagram, whic h can be explained by CH E. Since CHE is believed to be the key process to 
make a LGRB. lBestenlehner et al.l ( 120111 ) suggest VFTS 682 as a possible LGRB progenitor. 
A key ingredient to maintaining a CH E is sufficiently fast stellar rotation to induce a complete 
chemical mixing ( lSchwarzschildlll958[ ) . However, the problem is that CHE might cease when 
the strong mass-loss from VFTS 682 carries away too much angular momentum before its 
death. 

Theoretically low metallicity is favorable for single stars to be GRB progenitors, mainly 
because low metal licity leads to low mass-loss rates, that sustain the fast rotation required 
by LGRB models (lYoon fc Langer II2005I : IWoosley fc Heger 1120061 ). Recent simulation s pro- 
posed a metallicity threshold oi Z < 0.004 for GRB production ( lYoon et al.ll2006[ ). Ob- 
servations find that local LGRBs associated with supernovae (SNe) have me tal-poor host 
galaxies with oxygen abu ndance of 12 + log(0/H) < 8.6, oi Z < (0.2 — Q .b)Z(:^ ( 



20061 : iModjaz et al.ll2008l ). but not as metal-poor as required by models. IStanek et al 



Stane^ 



c et al. 



(120061 ) 
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argued that LGRBs trace only low-metallicit y star formation. Al though one or two high- 
metallicity GRB hosts have been discovered (ILevesque et al.ll2010l ). all of these GRBs lack 
an SN signature, leaving it unclear whether those GRBs can be linked with star formation 
or not. The half-solar metallicity of VFTS 682 is just in the upper metal range given by 
LGRB-SN events, so it is unclear whether it can be a LGRB progenitor just based on the 
metallicity criterion. A more direct constraint on the fate of VFTS 682 can be derived from 
its other physical parameters as follows. 

The observed present mass-loss rate of VFTS 6 82 is log(M/Mf;^ yf"''^^ ~ —4.4 ± 0.2, 
with an estimate present-day mass as M ~ ISOM© (IBestenlehner et al.ll201ll ). The mass- 
loss timescale for VFTS 682 of Twmd = MjM ~ 3.8 Myr is comparable to the stellar nuclear 
timescale Xnuc ~ 5 x 10^ yr {M /Mq){L/ Lq)~^ ~ 2.6 Myr for the current hydrogen abun- 
dance Xh = 0.55. This means that VFTS 682 can easily be stripped given its current 
mass-loss rate. For simplicity, if we take a remaining lifetime for VFTS 682 on the main- 
sequence (MS) as ~ 2 Myr with a constant mass-loss rate as its present value, the total 
mass lost when it evolves off the MS will be Mioss ~ 8OM0. Since the hydrogen abundance 
will continue to decrease fro m its current value, w hile the luminosity L changes a little with 
both deceasing M and Xh (iGrafener et all 1201 II ), the theoretically predicted mass- loss rate 
M oc 1O"° '^^^«(L/L0)° '^2 increases in the future jGrafener fc Hamannlbood ). If the hydro- 
gen abundance smoothly drops from Xh = 0.55 to Xh — 0, the total mass-loss increases 
from the rough estimate ~ 8OM0 to as high as ~ lOOM©. This would leave a remaining 
core Mr ~ 70 — 50Mq star as the star leaves the MS, and its subsequent evolution could 
eject even more mass. Unlike normal stars, which include core and envelope components, 
stars undergoing CHE ca n be approximated as a chemically mixed, rigidly rotating bodies 
( iMeynet fc Maederl 120001 ). Treating VFTS 682 as a rigid object with a current surface rota- 
tion velocity Vq and radius-mass relation R oc M", the final rotation velocity Vf at the end o f 
MS isVf/vQ = (Mr/Mo)*-^"^"^/^. Taking a ~ 1 for very massive stars (lYungelson et al.ll2008l ). 
we have Vf/vQ ~ 0.6 — 0.7 and the star will have lost more than 85% of its current angular 
momentum before its death. For currently vo < 700 km s~^, the stellar final rotation velocity 
Vf <500 km s~^ is unlik ely high enou gh to maintain C HE in the LMC (see the criterion in 



Mevnet fc Maederl l2000l . or Fig. 7 in iBrott et all 1201 If ). VFTS 682 will have a hydrogenic 



envelope at death, thus cannot become a LGRB. On the other hand, if VFTS 682 has an 
incredibly rapid rotation velocity vo > 700 km s ~^, the possibility o f maintaining CHE and 
producing a GRB associated with a hypernova ( iNomoto et al.ll2005l ) is not excluded. 

Some other effects such as shorter lifetimes or wind anisot ropies can decrease the a ngular 
momentum loss and help maintain the angular momentum ( iMeynet fc Maederl 120071 ). But 
even for a large helium star with sufficient angular momentum, there is still an extra problem 
for producing a GRB. A relativistic jet generated in the stellar center cannot break out of 
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the s tar if the duration of c entral engine is shorter than the jet crossing time inside the 
star (IMeszaros fc Reesll200ll ). Taking a radius of ~ 10^^ cm for a hehum star with a mass 
ry lOOMc:^, the jet cross ing time inside the star is estimated as tcmss ^ 100(rHp/10^^ cm) s 
(IMeszaros fc Reed 1200 ll ). which requires the duration of GRB central engine to be longer 
than 100 seconds. Such durat ions are observed, but rare, compared to those of durations 
< 100 s JSakamoto et aliboilh . 



4. Fate of Dusty Clouds and Possible Dark Gamma-Ray Burst? 

What is the fate of the dusty clouds near VFTS 682 after its death? Will the clouds 
be destroyed by the explosion of VFTS 682? Because the properties of dusty clouds around 
VFTS 682 are not at all un usual in the 30 Dor r egion, we adopt the gas density of the 30 Dor 
core region n ~ 20 cm~^ (IKawada et al.ll201ll ) with a typical value of Nn/Ay ~ 0.7 x 10^^ 



cm ^ in the LMC (jSchady et al.l 120071 ). This implies that the foreground extinction region 



has a size of ~ 50 pc for Ay ~ 4.45. This is consistent with the 8.0 fim image in Figure 
[H where the dusty region around VFTS 682 has a projected size of ~ 30 — 40 pc. In the 
"walkaway" scenario, VFTS 682 has both a tangential and RV velocity of ~30 km s~^ away 
from R136. Is so, VFTS 682 will probably escape the dusty region in its remaining ~ 2 Myr 
lifetime. On the other hand, if VFTS 682 formed in situ, or has a much shorter lifetime, it 
will be still in the very dusty clouds at its death. 

Since the luminosity of the afterglow of a GRB cannot be well determined from just 
the progenitor mass and rotation, it is uncertain whether the potential GRB from VFTS 
682 will be dark, or the dusty clouds will be destroyed by the optical flashes produced 
by VFTS 682. Typically we would assume a relativistic jet breaks out of the star and 
emits an isotropic gamma-ray luminosity Li^° ~ 10^^ erg s~^ (iLee et al.ll2000l ). followed by 

O.lLi' 



an early optical afterglow with a luminosity L 



ISO 

opt 



' ISO 

"7 • 



In this case any dust region 



smaller than R ~ 30{U^Z/10^^^ erg s ^y/"^ pc cannot survive because of dust sublimation 



by the optical flash an d fragmentation by the burst and afterglow (IWaxman fc Draind 12000 



Rei chart fc Pricdl2002l ) . In other words, the dusty region covering VFTS 682 and R136 would 
probably be destroyed due to an early optical afterglow L^p^. > 10^° erg s~^. Otherwise, the 
less luminous early X-ray/optical afterglows will heat and ionize the surrounding environment 
out to ~ 100 pc, decreasing the dust column density on a timescale of tens to hundreds of 
minutes in the observer's frame (jPerna fc Loebl 119981 ). 



However, it will be another story if VFTS 682 produces an intrinsically faint intrinsic 
faint GRB. It is possible that the inefficient jet propagation inside the star only gives an 
underluminous burst, as well as a dim optical afterglow {L'-oZ < 10^° erg s~^) linked with the 
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weak burst. Then the dusty clouds will be heated and ionized, but not be totally destroyed. 



Conclusions and Discussion 



We give a first order estimate whether VFTS 682 can be a LGRB progenitor in the sense 
of CHE, which is still a question mark in previous work. VFTS 682 will most likely lose more 
than 85% of its current angular momentum and cease its CHE evolution before it leaves MS 
with a remaining mass ~ 50 — 70 Mq. In general VFTS 682 will fail to produce a GRB due 
to its strong mass-loss and the resulting angular momentum loss, unless it has a currently 
extreme rapid rotation [vq > 700 km s~^) which helps it maintain CHE and eventually 
produce a GRB associated with a hypernova. Wind anisotropics and shorter lifetime could 
leave a more massive faster rotating star, but it is doubtful whether a relativistic jet could 
travel through such a thick stellar envelope and break out of the stellar surface, unless it is 
a rare long-lived (> 100 s) LGRBs. Similarly, it is unlikely that VFTS 682 will be heavily 
obscured at death and produce a dark GRB. Its proper motion will probably cause it out of 
the dusty region, and a GRB of an early optical afterglow Lopt > 10^° erg s~^ would destroy 
the dust clouds within 30 pc, otherwise the dusty clouds can be heated and ionized up to a 
region of ~ 100 pc by X-ray and optical radiation given by the death of VFTS 682. 

CHE is believed to be the crucial process in the evolution path towards LGRBs. How- 
ever, the observed sample of WRs likely undergoing CHE is quite small. The only observation 
of WR stars besides VFTS 68 2 which might be un dergoing are two WNh stars in the SMC 
(i.e., SMC-WRl and WR2 in Martins et al.ll2009l ). which makes these three WR stars ex- 
tremely important to understanding the evolution of WR stars undergoing CHE and the 
related problem of LGRB formation. The angular momentum losses of the two WR stars 
in the SMC will be much le s s sign ificant than VFTS 682 in the LMC. Taking the stellar 
parameters in iMartins et al.l (120091 ) . we estimate that the nuclear timescale of SMC-WRl 
(WR2) Tnuc ~ 4 Myr (5 Myr) is shorter than the wind timescale r^md ~ 10 — 20 Myr, and 
the final stellar rotation velocity should be 90% of the current velocity. Therefore, the two 
WR stars in the SMC are more likely GRB progenitors in t he scenario of CHE, although the 
metallicity threshold is still an issue ( iMartins et al.ll2009l ). In any case, finding CHE WR 
stars in the LMC and SMC promotes a future work on theory models. 

VFTS 682 is a very interesting star. As mentioned in Section 2, because there is no 
observation evidence to show that VFTS 682 is unusual compared to other stars in the 30 
Dor region, the interesting possibility of existence of others massive stars like VFTS 682 
in the vicinity of R136 is not excluded. Note that R136 is sufficie nt young and massive 
(< 5.5 X IO^Mq) to generate runaway stars beyond 15OM0. Recently iBanerjee et al.l (1201 ll ) 
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gives a theoretical model study on the dynamical ejection of runaway massive stars from 
R136. We suggest that there might be other massive stars that "walked away" from R136, 
but are currently hidden behind even more dust than VFTS 682. Since mid-IR date cannot 
be used to flag such stars, spectroscopic observations, which are beyond the scope of this 
Letter, should be further investigated to show the possibility of their existence. 

We are grateful to C. S. Kochanek, T. A. Thompson and S. Horiuchi for useful comments 
and critical reading of the manuscript, and M. Pinsonneault, R. Khan, J. Beacom and J. 
Skowron for helpful discussions. This work is partially supported by the NSF grant AST- 
1108687. 
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